Astronomy & Astrophysics manuscript no. w28 aa 


© ESO 20081 


February 2, 2008 





Discovery of very high energy gamma-ray emission coincident with 
molecular clouds in the W 28 (G6.4-0.1) field 

F. Aharoniani'i^, A.G. Akhperjanian^, A.R. Bazer-Bachi^', B. Behera'"^, M. Beilicke^, W. Benbow', D. Berge^ *, 
K. Bernlohr''^ C. Boisson^, O. Bolz\ V. Borrel^ I. Braun^ E. Brion'', A.M. Brown^ R. Buhler^ T. Bulik^^, 
I. Busching'^, T. Boutelier^^, S. Carrigan^, P.M. Chadwick^ L.-M. Chouneti", A.C. Clapson', G. Coignet'i, 
R. Cornils'*, L. Costamante^-^^, B. Degrange'", H.J. Dickinson^, A. Djannati-Atai'^^, W. Domainko\ L.O'C. Drury^^, 
G. Dubus'°, J. Dyks^"*, K. Egberts', D. Emmanoulopoulos'^, R Espigat'^, C. Farnier'^, F. Feinstein'^, A. Fiasson'^, 

A. Forster', G. Fontaine'°, Y. Fukui^^, Seb. Funk^ S. Funk', M. FuBling^ Y.A. Gallant'^ B. Giebels'°, 
J.F Glicenstein'^, B. Gluck'^, R Goret'', C. Hadjichristidis^, D. Hauser', M. Hauser''^, G. Heinzelmann"', G. Henri'"', 
G. Hermann', J.A. Hinton'-"^ **, A. Hoffmann'^ W. Hofmann', M. Holleran'^, S. Hoppe', D. Horns'^ 
A. Jacholkowska'^, O.C. de Jager^, E. Kendziorra'*^, M. Kerschhaggl^ , B. Khelifi'"'', Nu. Komin'^, K. Kosack', 
G. Lamanna", I.J. Latham^, R. Le Gallou^, A. Lemiere'^, M. Lemoine-Goumard'", J.-R Lenain^, T. Lohse^, 
J.M. Martin^, O. Martineau-Huynh'^, A. Marcowith'^, C. Masterson'^, G. Maurin'^, T.J.L. McComb^ R. Moderski^"^, 
Y. Moriguchi26, E. Moulin'^ ^ M. de Naurois'^ D. NedbaP, S.J. Nolan^ J-R Olivet K.J. Orford^ J.L. Osborne^ 
M. Ostrowski^^, M. Panter', G. Pedaletti'", G. Pelletier'^, R-0. Petrucci''', S. Pita'^, G. Puhlhofer'^, M. Punch'^, 
S. Ranchon", B.C. Raubenheimer*^ , M. Raue", S.M. Rayner^ O. Reimer M. Renaud', J. Ripken"*, L. Rob^o, 
L. Rolland^, S. Rosier-Lees", G. Rowell' \ B. Rudak^", J. Ruppel^', V. Sahakian^, A. Santangelo'^ L. Sauge'^, 
S. Schlenker^, R. Schlickeiser^' , R. Schroder^', U. Schwanke^, S. Schwarzburg'^, S. Schwemmer'^, A. Shalchi^', 
H. Sol^, D. Spangler*^, L. Stawarz^^, R. Steenkamp^^, C. Stegmann'^, G. Superina'°, T. Takeuchi^^, PH. Tam"', 
J.-P Tavernet''^, R. Terrier'^, C. van Eldik', G. Vasileiadis'^ C. Venter'^, J.P Vialle", P Vincent'^, M. Vivier^, 

H.J. Volk', F Volpe'", S.J. Wagner'", and M. Ward^ 

Correspondence to:growell@physics.adelaide.edu.au, mori@aserv.a.phys.nagoya-u.ac.jp 

(Affiliations can be found after the references) 
Received / Accepted 

ABSTRACT 

Aims. Observations of shell-type supernova remnants (SNRs) in thie GeV to multi-TeV y-ray band, coupled with those at millimetre radio 
wavelengths, are motivated by the search for cosmic-ray accelerators in our Galaxy. The old-age mixed-morphology SNR W 28 (distance ~2 kpc) 
is a prime target due to its interaction with molecular clouds along its northeastern boundary and other clouds situated nearby. 
Methods. We observed the W 28 field (for ~40 h) at very high energy (VHE) y-ray energies (E >0.1 TeV) with the H.E.S.S. Cherenkov telescopes. 
A reanalysis of EGRET E >100 MeV data was also undertaken. Results from the NANTEN 4m telescope Galactic plane survey and other CO 
observations were used to study molecular clouds. 

Results. We have discovered VHE y-ray emission (HESS J1801-233) coincident with the northeastern boundary of W 28 and a complex 
of sources (HESS J1800-240A, B and C) -0.5° south of W 28 in the Galactic disc. The EGRET source (GRO J1801-2320) is centred on 
HESS J1801-233 but may also be related to HESS J1800-240 given the large EGRET point spread function. The VHE differential photon spectra 
are well fit by pure power laws with indices T ~ 2.3 to 2.7. The spectral indices of HESS J1800-240A, B, and C are consistent within statistical 
errors. All VHE sources are ~10' in intrinsic radius except for HESS J1800-240C, which appears pointlike. The NANTEN '^CO(J=1-0) data 
reveal molecular clouds positionally associating with the VHE emission, spanning a ~15 km s"' range in local standard of rest velocity. 
Conclusions. The VHE/molecular cloud association could indicate a hadronic origin for HESS J1801-233 and HESS J1800-240, and several 
cloud components in projection may contribute to the VHE emission. The clouds have components covering a broad velocity range encompassing 
the distance estimates for W 28 (~2 kpc) and extending up to ~4 kpc. Assuming hadronic origin and distances of 2 and 4 kpc for cloud components, 
the required cosmic-ray density enhancement factors (with respect to the solar value) are in the range ~10 to ~30. If situated at 2 kpc distance, 
such cosmic-ray densities may be supplied by SNRs like W 28. Additionally and/or alternatively, particle acceleration may come from several 
catalogued SNRs and SNR candidates, the energetic ultra compact HII region W 28A2, and the HII regions M 8 and M 20, along with their 
associated open clusters. Further sub-mm observations would be recommended to probe in detail the dynamics of the molecular clouds at velocites 
>10 km s"' and their possible connection to W 28. 
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1. Introduction: W 28 and surroundings 

The study of shell-type supernova remnants (SNRs) at y-ray en- 
ergies is motivated by the long-held idea that they are the dom- 
inant sites of hadronic Galactic cosmic-ray (CR) acceleration 
to energies approaching the knee (~ lO'^ eV) (e.g. Ginzburg 
& Syrovatskii 1964, Blandford & Eichler 1987). CRs (hadrons 
and electrons) are injected into the SNR shock front, and are 
then accelerated via the diffusive shock acceleration (DSA) pro- 
cess (for a review see Drury 1983). Subsequent y-ray production 
from the interaction of these CRs with ambient matter and/or 
electromagnetic fields is a tracer of such non-thermal particle 
acceleration, and establishing the hadronic or electronic nature 
of the parent CRs in any y-ray source remains a key issue. Two 
SNRs, RX J1713.7-3946 and RX J0852.0-4622, have so far es- 
tablished shell-like morphology in VHE y-rays (Aharonian et al. 
2004a, 2005c, 2006b, 2007b, 2007a), with spectra extending to 
20 TeV and beyond. In particular for RX J1713.7-3946, par- 
ticle acceleration up to at least 100 TeV is inferred from the 
H.E.S.S. observations. Although a hadronic origin of the VHE y- 
ray emission is highly likely in the above cases (Aharonian et al. 
2006b, 2007b, Berezhko & Volk 2006, Berezhko, Puhlhofer & 
Volk 2007), an electronic origin is not ruled out. 

Disentangling the electronic and hadronic components in 
TeV SNRs may be made easier by studying: (1) SNR y-ray spec- 
tra well beyond ~ 10 TeV, an energy regime where electrons suf- 
fer strong radiative energy losses and due to Klein-Nishina ef- 
fects the resulting inverse-Compton spectra tend to show a cut- 
off; (2) older SNRs (age approaching 10^ yr) in which accel- 
erated electrons have lost much of their energy through radia- 
tive coohng and do not reach multi-TeV energies; (3) SNRs in- 
teracting with adjacent molecular clouds of very high densities 
n > 10' cm"'. It is the latter regard especially (and to a certain 
degree the second) which makes the SNR W 28 (G6.4-0.1) an 
attractive target for VHE y-ray studies. In this paper we outline 
the discovery of VHE y-ray emission from several sites in the 
W 28 field and briefly discuss their relationship with molecu- 
lar clouds, W 28, and other potential particle accelerators in the 
region. 

W 28 (G6.4-0. 1) is a mixed-morphology SNR, with dimen- 
sions 50'x45' and an estimated distance between 1.8 and 3.3 kpc 
(eg. Goudis 1976, Lozinskaya 1981). It is an old-age SNR (age 
35000 to 150000 yr; eg. Kaspi et al. 1993), thought to have en- 
tered its radiative phase of evolution (eg. Lozinskaya 1981) in 
which much of its CRs have escaped into the surrounding inter- 
stellar medium (ISM). We note also that the evolutionary status 
(Sedov and/or radiative) of shell-type SNRs may depend on the 
density of their surroundings (see eg. Blondin et al. 1998). 

W 28 is distinguished by its interaction with a molecular 
cloud (Wootten 1981) along its north and northeastern bound- 
aries. This interaction is traced by the high concentration of 
1720 MHz OH masers (Frail et al. 1994, Claussen et al. 1997, 
1999), and also the location of very high-density (n > lO"' cm"') 
shocked gas (Arikawa et al. 1999, Reach et al. 2005). The shell- 
Uke radio emission (Long et al. 1991, Dubner et al. 2000) peaks 
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at the northern and northeastern boundaries where interaction 
with the molecular cloud is estabhshed. Further indication of the 
influence of W 28 on its surroundings is the expanding HI void 
at a distance ~1.9 kpc (Velazquez et al. 2002). The X-ray emis- 
sion, which overall is well-explained by a thermal model, peaks 
in the SNR centre but has local enhancements in a region over- 
lapping the northeastern SNR/molecular cloud interaction (Rho 
& Borkowski 2002). 

In the neighbourhood of W 28 are the radio-bright HII re- 
gions M 20 (Trifid Nebula at d -1.7 kpc Lynds etal. 1985 
- with open cluster NGC 6514), M 8 (Lagoon Nebula at 
d ~ 2 kpc TothiU etal. 2002 — containing the open clus- 
ters NGC 6523 and NGC 6530) and the ultra-compact HH 
region W 28A2, all of which are representative of the mas- 
sive star formation taking place in the region. Further discus- 
sion concerning the active star formation in this region may be 
found in van den Ancker et al. (1997) and references therein. 
Additional SNRs in the vicinity of W 28 have also been iden- 
tified: G6.67-0.42 and G7.06-0.12 (Yusef-Zadeh et al. 2000), 
G5.55-i-0.32, G6.10-h0.53 and G7.20-i-0.20 (Brogan et al. 2006). 
The pulsar PSR J1801-23 with spin-down luminosity E ~ 
6.2 X lO'"* erg s"' and distance d - 13.5 kpc (based on its disper- 
sion measure) is at the northern radio edge (Kaspi 1993). More 
recent discussion (Claussen et al. 2002) assigns a lower limit of 
9.4±2.4 kpc for the pulsar distance. 

W 28 has also been linked to y-ray emission detected at 
E > 300 MeV by COS-B (Pollock 1985) and E > 100 MeV by 
EGRET (Sturner & Dermer 1995, Esposito et al. 1996, Zhang 
etal. 1998). The EGRET source, listed in the 3rd catalogue 
(Hartman et al. 1999) as 3EG J1800-2338, is positioned at flie 
southern edge of the radio shell. We have also performed an anal- 
ysis of EGRET data, with additional data not included in the 3rd 
catalogue, and results are discussed later in this paper. 

Previous observations of the W 28 region at VHE ener- 
gies by the CANGAROO-I telescope revealed no evidence for 
such emission (Rowell et al. 2000) and upper limits at the 
~0.2 to 0.5 Crab-flux level for energies E > 1.5 TeV (1.1 to 
2.9x10"^' erg cm"^ s ') were set for various regions. 

2. Results at VHE and E >100 MeV 7-ray energies 

2.1. H.E.S.S. VHE analysis and results 

The High Energy Stereoscopic System (H.E.S.S.) was used to 
observe the W 28 region. Operating in the Southern Hemisphere, 
H.E.S.S. consists of four identical 13 m diameter Cherenkov 
telescopes (Bemlohr et al. 2003). H.E.S.S. employs the stereo- 
scopic imaging atmospheric Cherenkov technique, and is sensi- 
tive to y-rays above an energy threshold of ~0.1 TeV (Funk et al. 
2004). An angular resolution of 5' to 6' (Gaussian standard devi- 
ation) on an event-by-event basis is achieved, and the large field 
of view (FoV) witii fuU width at half maximum FWHM ~ 3.5° 
permits survey coverage in a single pointing. A point source 
sensitivity approaching 0.01 Crab flux (~ 10"'-' erg cm^^s"' 
at 1 TeV) is achieved for a 5cr detection after ~25 hr observa- 
tion. Further details concerning H.E.S.S. can be found in Hinton 
(2004) and references therein. 

The total observation time covering the W 28 region amounts 
to ~42 hr in a series of runs (with typical duration ~28 min) 
spread over the 2004, 2005 and 2006 seasons. Runs were ac- 
cepted for analysis if they met quality control criteria based on 
the recorded rate of isotropic CR background events, the number 
of malfunctioning pixels in each camera, the calibration and the 
tracking performance (see Aharonian et al. 2004b for details). 
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Data were analysed using the moment-based Hillas anal- 
ysis procedure, the same used in the analysis of the inner 
Galactic Plane Scan datasets (Aharonian et al. 2005a, 2006a). 
Observations covered a range of zenith angles leading to energy 
thresholds of ~ 320 Ge V with hard cuts (Cherenkov image inte- 
grated intensity or size >200 photoelectrons) and ~150 GeV for 
standard cuts (size > 80 photoelectrons). Hard cuts were used in 
VHE y-ray images, source location studies and energy spectra. 
In addition. Standard cuts were used in energy spectra in order 
to increase the energy coverage of extracted spectra. Generally 
consistent results were obtained using an alternative analysis 
based on a model of Cherenkov image parameters (de Naurois 
2006), which also utilises an independent calibration and lower 
cut on image size of >60 photoelectrons. A forthcoming paper 
will highhght results in detail from this analysis, which achieves 
improved sensitivities at lower thresholds compared to the pure 
Hillas-based analysis. 

The VHE y-ray image (Fig. 1) reveals two sites of VHE y- 
ray emission in the direction of the northeastern and southern 
boundaries of the W 28 SNR. The colour scale in this figure 
depicts the Gaussian-smoothed VHE excess counts above a CR 
background estimate according to the template model (Rowell 
2003), along with significance contours obtained after integrat- 
ing events within a radius of 0.1° from each bin centre (appro- 
priate for pointlike source searching). Similar images were ob- 
tained using alternative CR background models. A smoothing 
radius of 4.2' was used to sufiiciently smooth out random fluc- 
tuations in the image. An assessment of the VHE post-trial sig- 
nificances was made from our original search for marginally ex- 
tended sources, which employed an a priori integration radius 
9 = 0.2°. Under this scheme we applied ~ 2.2 x 10^ trials (a very 
conservative value applied to these data) accumulated in search- 
ing for sources in the inner Galactic Plane (as in Aharonian 
2005a). The pre-trial significance of the VHE sources, at > +la, 
is therefore converted to a post-trial significance of > +5cr. 

Based on the significance contours in Fig. 1, we assign la- 
bels to the northeastern source, HESS J1801-233, and to the 
complex of sources to the south, HESS J 1800-240, according 
to their best fit positions (fitting a 2D Gaussian and ellipse re- 
spectively to the unsmoothed excess map). Three components 
of HESS J1800-240 are identified, labeled here A, B and C 
from East to West. These components represent local peaks 
~ 2cr above their surrounds. Although not convincingly resolved 
under this analysis these components may comprise separate 
sources (or at least in part) due to their possible relationship with 
distinct multiwavelength counterparts (discussed later). 

Difl'erential photon energy spectra were extracted 
from HESS J1801-233 and all three components of 
HESS J1800-240. Spectra were well-fit by pure power 
laws {cLNIdE = kiE/lTeV)-^) with photon indices F ~ 2.5 
to 2.7 in the energy range ~0.3 to ~5 TeV (see Table 1 for 
results). Spectral fits were obtained using fluxes from a com- 
bination of hard and standard cuts to maximise the energy 
coverage. Spectral analysis employed the reflected background 
model (Berge etal. 2007), in which control regions reflected 
through each tracking position (taking care to avoid known 
VHE y-ray sources) were used to estimate the CR background. 
Within the statistical and systematic errors, the photon indices 
appear consistent throughout HESS J1800-240. Except for 
HESS J1800-240C, aU of flie VHE sources appear extended 
with intrinsic radii of ~10'. At a distance of 2 kpc, the VHE 
source luminosities in the energy range 0.3 to 3 TeV would be 
on the order of 10^^ erg s"^ . 



2.2. EGRET E > 100 MeV analysis and results 

We have also analysed EGRET data for the W 28 region, using 
CGRO observation cycles (OC) 1 to 6. This shghtly expands 
on the dataset of the 3rd EGRET catalogue (using OCs 1 to 
4; Hartman et al. 1999), which revealed the pointlike source, 
3EG J1800-2338 (E > 100 MeV). Our analysis confirms 
the presence of a pointlike £ > 100 MeV source in this re- 
gion, here labeled GRO J1801-2320 (for E > 100 MeV). 
GRO J1801-2320 appears slightly shifted (-0.2°) with respect 
to the 3EG position. The 3EG position refers to a £ >100 
MeV determination based on the difiiise model as of Hunter 
etal. (1997). Our dedicated analysis of archival EGRET data 
comprises different analysis compared to the 3EG catalogue. 
We first employed the finalised EGRET instrumental responses, 
which were made available by 2001 and are considered manda- 
tory for investigating an EGRET source under conditions ap- 
plicable from the end of OC 4 (narrow field of view modus; 
rapidly deteriorating spark chamber eflficiency; and other issues). 
Second, we restricted the analysis both in narrowing the data 
selection to pointing angles with respect to our region of inter- 
est, which avoids the need to invoke a wide-angle point spread 
function (PSF). Thirdly, the imprecision of the interstellar emis- 
sion model was countered via adjustments on analysis parame- 
ters gmult and gbias to account for local deviations from the 
large-scale diffuse emission model in the region of interest. The 
68% and 95% location contours of GRO J1801-2320 are plot- 
ted in Fig. 1, and match weU the location of HESS J1801-233. 
Since however the EGRET degree-scale PSF easily encom- 
passes both of the VHE sources, we cannot rule out a rela- 
tionship with HESS J1800-240. For the energy spectrum of 
GRO J 180 1-2320, we have used the flux points extracted at 
the position of 3EG J1800-2338 as negligible differences were 
found between ours and that obtained at the nominal 3EG po- 
sition. Fitting a pure power law we obtained a spectral index of 
F = 2.16 + 0.10, quite consistent with the published value from 
Hartman etal. (1999). Comparisons are made with the VHE 
spectrum of HESS J1801-233 and HESS J1800-240 in §5. 

3. NANTEN and other observations of Molecular 
Clouds 

In searching for molecular cloud counterparts to the VHE 
sources, we analysed '^CO (/=l-0) molecular line observa- 
tions taken by the 4-meter mm/sub-mm NANTEN telescope, 
at Las Campanas Observatory, Chile (Mizuno & Fukui 2004). 
The NANTEN Galactic Plane Survey data of 1999 to 2003 (see 
Matsunaga et al. (2001) and references therein for details) were 
used, and for the W 28 region, the survey grid spacing was 4'. 

Figure 2 (upper left panel) shows the '^CO (/=l-0) image 
integrated over the Local Standard of Rest velocity (Vlsr) range 
to 10 km s"' , while the right panel shows the image integrated 
over the range Vlsr=10 to 20 km s"^ Two prominent '^CO 
features representing molecular clouds centred at (l,b)=(6.7°, 
-0.3°) and (/,fo)=(5.9°, -0.4°) spatiafly correspond with the 
VHE y-ray emission. As shown in Fig. 2, these molecular clouds 
span both Vlsr ranges. According to the Galactic rotation model 
of Brand & Blitz (1993), these Vlsr ranges formally correspond 
to kinematic distances of approximately to ~2.5 kpc (over- 
lapping the Sagittarius arm), and 2.5 to ~4 kpc (reaching the 
Scutum-Crux arm) respectively. Given the uncertainties in rota- 
tion models close to the Galactic centre, such Vlsr ranges would 
cover the distance estimates for W 28, the most prominent SNR 
in the region. Much discussion has centred on the systemic ve- 
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Fig.l. Image (1.5° x 1.5°) of 
the VHE y-ray excess counts 
(events), corrected for exposure 
and smoothed with a Gaussian of 
radius 4.2' (standard deviation). 
Overlaid are solid green contours 
of VHE excess (pre-trial) signif- 
icance levels of 4, 5, and 6cr, af- 
ter integrating events within an 
oversampling radius 0=0.1° ap- 
propriate for pointlike sources. 
The thin-dashed circle depicts 
the approximate radio boundary 
of the SNR W 28 guided pre- 
dominantly by the bright north- 
ern emission (see Dubner et al. 
2000 & Brogan et al. 2006). 
Identified here are VHE source 
regions HESS J1801-233 to 
the northeast, and a complex 
of sources HESS J 1 800-240 
(A, B & C) to the south of 
W 28. Also indicated are: HII 
regions (black stars); W 28 A2 
(see text), G6. 1-0.6 (Kuchar 
& Clark 1997), 6.225-0.569 
(Lockman 1989); The 68% and 
95% location contours (thick- 
dashed yellow lines) of the 
£ > 100 MeV EGRET source 
GRO J1801-2320; the pulsar 
PSR J1801-23 (white triangle). 
The inset to the bottom left 
depicts a pointlike source for 
this analysis after the Gaussian 
smoothing applied to the main 
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Table 1. Numerical summary^ * for the VHE and E > 100 MeV sources in the W 28 region including positional and spectral information. 



Best fit position (J2000.0) Spectral analysis 



Name 




R.A. [deg] 




Dec [deg] 




'cTsrc 


[deg] 


[a-] (evts) 








4p 






HESS J1801 


-233 


270.426 ± 


0.031 


-23.335 ± 


0.032 


0.17 


±0.03 


-H7.9 (281) 


7.50 ± 


1.11 ± 


0.30 


2.66 


± 0.27 


1.5 


HESS J 1800 


-240A§ 


270.491 ± 


0.001 


-23.962 ± 


0.001 


0.15 




-H6.0 (180) 


7.65 ± 


1.01 ± 


0.50 


2.55 


±0.18 


1.5 


HESS J 1800 


-240B§ 


270.110 ± 


0.002 


-24.039 ± 


0.009 


0.15 




-1-7.8 (236) 


7.58 ± 


0.90 ± 


0.15 


2.50 


± 0.17 


1.4 


HESS JI800 


-240C 


269.715 ± 


0.014 


-24.052 ± 


0.006 


0.02 


±0.15 


+4.5 (71) 


4.59 ± 


0.89 ± 


0.20 


2.31 


±0.35 


0.8 


HESS JI800 


-240^5 


270.156 ± 


0.044 


-23.996 ± 


0.022 


0.321^ ± 0.05 


+ 10.3 (652) 


18.63 


± 1.85 


t 1.20 


2.49 


± 0.14 


3.6 














0.17 


0='^ ± 0.03 
















GRO JI801- 


2320 


270.360 ± 


0.150 


-23.340 ± 


0.150 






+ 13.2 


3.35 ± 


0.52 




2.16 


±0.10 


480.0 



t VHE photon spectra are derived from a region of radius 6 - ^|0A^~^o^ centered on each source's position unless otherwise indicated. 
t In spectra, a function dN/dE = kE^^ ph cm"" s"' TeV"' is fitted. E is in TeV units (H.E.S.S. data); GeV units (EGRET data). 

1. Fitted intrinsic source size (Gaussian std. dev.) 

2. Statistical significance and excess events in brackets; for H.E.S.S. sources using Li & Ma (1983); for EGRET sources given by 5 - ^/J\ for defined by Mattox et al. (1996) 

3. For H.E.S.S. sources: xlO"'^ ph cm"" s"^ TeV"' at 1 TeV (with statistical and systematic errors); For EGRET sources; xlO"^ ph cm"" s"' GeV"' at 1 GeV (with statistical errors) 

4. Only statistical errors indicated. Systematic eiror is estimated at ±0.2 

5. Luminosity X 10'^ erg s"' at 2 kpc (0.3 to 3 TeV for H.E.S.S.; 0.04 to 6 GeV for EGRET) 

§Due to cross contamination between components A & B, a fixed value of tr^rc - 0.15° estimated visually from Fig. 1 was used. 

§§Spectrum extracted from a 0.8° X 0.6° elliptical region encompassing all components A, B, C, and matching the size of the coiresponding molecular cloud. 



locity (SV) of W 28 (and hence its distance), and how much 
W 28 has influenced matter in the region. Ha (Radhakrishman 
et al. 1972) and HI absorption features (Lozinskaya et al. 1981) 
have suggested SV~18 km s"'. Claussen et al. (1997) have 
pointed to SV~17 km More recent HI studies by Velazquez 
et al. (2002) suggest SV=-h7 km s"' (which leads to the distance 
estimate for W 28 at ~ 1.9 kpc). They also suggest a HI shell 



may also extend over the Vlsr = -25 to -t-38 km s"' range, giv- 
ing rise to a shock speed of ~30 km s"' . Torres et al. (2003) and 
Reach et al. (2005) have also studied the large-scale '^CO(J=l- 
0) emission for this region using the survey data of Dame et al. 
(2001), suggesting that the parent molecular cloud under the 
influence of W 28 is presently centred at Vlsr ~19 km s"'. 
The Galactic longitude-velocity (I-v) diagram (bottom panels of 
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Fig. 2) from our NANTEN data integrated over the Galactic lat- 
itude ranges b = -0.125° to -0.5° andb = -0.125° to -0.7° 
shows the distribution of molecular material in relation to the 
SV of W 28 from the HI studies of Velazquez. The wider, lat- 
ter b range shows the effect of including the cloud component 
overlapping HESS J1800-240A. A void or dip in CO emission 
appears at a similar Vlsr range as found in the HI data, with 
much of the molecular material appearing to surround the void 
in positive Vlsr values with respect to the SV of W 28. A similar 
longitude- velocity picture was revealed by Torres et al. (2003) 
(see their Fig. 22). 

The Vlsr=0 to 10 km s~^ component of the northeast 
cloud overlapping HESS J1801-233 is already well studied (see 
Reach et al. 2005 and references therein). Shocked '^CO(7=3- 
2) molecular gas as indicated by a broad wing-like line disper- 
sion (Arikawa et al. 1999 — hereafter A99; using the James 
Clerk Maxwell Telescope (JCMT); in 15" grid steps) and a high 
concentration of OH masers (Claussen et al. 1997), suggests ma- 
terial here has been compressed by the SNR shock in W 28. The 
line dispersion, AV < 70 km s"\ is an indicator of the SNR 
shock speed in this particular region. The unshocked gas was 
also mapped by A99 via '^CO(7=1-0) observations with the 
Nobeyama 45 m telescope (in 34" grid steps for Vlsr - +4 
to +9 km s~^). The shocked and unshocked gas extends to the 
northeast and northern boundaries of W 28 (see Fig. 3 of A99), 
and it appears just their northeastern components are position- 
ally coincident with the VHE emission of HESS J1801-233. 
A99 estimate the mass and average density of the shocked gas 
at M ~ 2 X 10^ Mq and n ~ 10"* cm"^ respectively. For the un- 
shocked gas, A99 obtained M ~ 4 x 10^ M© and n ~ 10^ cm"^ 
respectively. The Vlsr=10 to 20 km s~' range in our NANTEN 
data also reveals additional molecular clouds along the line of 
sight that could contribute to the VHE emission. 

The southern cloud overlaps all components of 
HESS J1800-240, with a dominant fraction of the cloud 
overlapping components A and B. The component of this cloud 
visible in the Vlsr = to 10 km s~^ range coincides well with 
HESS J1800-240B and the HII region W 28A2. The strongest 
CO temperature peak of this component at (/, Z?)=(5.9°, -0.4°) 
is within 0.02° of W 28A2, and is likely the dense material 
surrounding this Hll region. Moreover the peak's velocity at 
Vlsr 9-10 km s"' (with dispersion of ~15 km s"'), suggests 
a distance (~2.4 kpc) similar to that of W 28A2 (~2 kpc; 
Acord et al. 1997), and also W 28. In the Vlsr = 10 km s'^ to 
20 km s"' range, molecular material appears to coincide with 
all three VHE components of HESS J 1800-240. In particular, 
HESS J1800-240A and C have molecular cloud overlaps only 
in this latter Vlsr range. 

Using the relation between the hydrogen column density 
A'(H2) and the '2cO(/=l-0) intensity (the X-factor) ^('^00), 
NiU2) = 1.5 X 102" [W(i2C0)/(Kkm/s)] (cm'^) (Strong etal. 
2004), we estimate a total mass for the northeastern cloud from 
our NANTEN data at ~ 5 x 10"* Mq for d =2 kpc within an 
elliptical region of diameter 0.2° x 0.4° (7x14 pc; centred on 
HESS J1801-233) for the velocity range 0-25 km s '. An av- 
erage density (for neutral hydrogen) of ~ 1.4x10^ cm"^ is also 
derived. Similarly the total mass of the southern cloud is es- 
timated at ~ 1.0x10^ Mq for d=2 kpc and combining clouds 
from a circular area of radius 0.15° (5 pc) for the velocity range 
12-20 km s \ and area 0.3° x 0.6° (10.5x21 pc) in diame- 
ter for the velocity range 0-12 km s"^ (both regions are cen- 
tred on HESS J1800-240B). The corresponding average density 
is ~1.0xl0-' cm"-'. By integrating over the rather broad 0-20 
km s ' and 0-25 km s"^ ranges we assume that the molecular 



material along this line of sight is physically connected at the 
same distance (for example d ~2 kpc) and possibly distrupted 
or shocked by a local energy source. Systematic effects in the 
mass estimates arise from the velocity crowding in this part of 
the Galactic plane, and also the broad velocity range for which 
X-factor used above may not necessary apply. In the latter case, 
the X-factor may underestimate the cloud mass since an appre- 
ciable fraction of gas may be heated under the assumption of 
distrupted and/or shock-heated gas. One must allow for ~4 kpc 
distances for some or even all of the Vlsr >10 km s"' cloud 
components, and therefore the conclusion that they are not re- 
lated to W 28 and other interesting objects at d ~ 2 kpc. If the 
clouds are related, W 28 could play a disrupting role. The level 
of this disruption is however unclear since several other plausible 
candidates related to the star formation (discussed later) in this 
region could also contribute. Some other molecular cloud com- 
plexes have also been discussed as possibly disrupted by adja- 
cent SNRs and/or energetic sources (eg. Yamaguchi et al. 1999, 
Moriguchi et al. 2000). In table 2, we present a full summary 
of cloud masses and densities (for regions centered on the VHE 
source coordinates as in table 1) for various combiniations of 
cloud components and distances of 2 and 4 kpc. Velocity separa- 
tion of cloud components are based on their apparent distribution 
in Fig. 2 (bottom panels). 

4. Radio to X-ray views 

Figure 3 compares the radio (left panel), infrared and X-ray 
views (right panel) of the W 28 region with the VHE significance 
contours. The Very Large Array (VLA) 90 cm continuum radio 
image from Brogan et al. (2006) illustrates the shell-like SNR 
morphology peaking strongly along the northern and eastern 
boundaries. HESS J1801-233 can be seen to overlap the north- 
eastern shell of the SNR, coinciding with a strong peak in the 
90 cm continuum emission. We note that a thermal component is 

likely present in this peak, given its spectral index a 0.2 (for 

S oc v") between 90 and 20 cm (Dubner et al. 2000). Outlines 
of the SNRs traced by non-thermal radio emission, G6.67-0.42 
and G7.06-0.12 (Yusef-Zadeh etal. 2000, Helfand etal. 2006, 
labelled as G6.5 1-0.48 and G7.0-0.1 by Brogan etal. 2006) 
are also indicated. In addition, Brogan et al. notes that the 
non-thermal radio arc G5.7 1-0.08, which overlaps well with 
HESS J1800-240C, could be a partial shell and therefore an 
SNR candidate. The distances to G6.67-0.42 and G5.71-0.08 
are presently unknown. Directly south of W 28, the ultracom- 
pact HII region W 28A2 is a prominent radio source, and is 
positioned within 0.1° of the centroid of HESS J1800-240B. 
The other HII regions G6.1-0.6 (Kuchar & Clark 1997) and 
6.225-0.569 (Lockman 1989) are also associated with radio 
emission. 

The X-ray morphology as shown (Fig. 3 right panel) in the 
ROSAT PSPC (0.5 to 2.4 keV) image from Rho & Borkowski 
(2002) reveals the central concentration of X-ray emission, 
which is predominantly thermal in nature with characteristic 
temperatures in the range kT -0.4 to 2 keV. An X-ray peak or 
Ear Ues at the northeastern boundary and just outside the 4cr 
significance contour of HESS J1801-233. A non-thermal com- 
ponent to the ear emission (3x1.5') (2.1x10"'"* erg cm"^ s"' at 
1 keV) with a power-law index r=1.3 has been suggested by 
Ueno et al. (2003a) based on XMM-Newton observations in the 
0.5 to 10 keV energy range. The total kinetic energy of the SNR 
is estimated at ~ 4 x 10^'' erg, which could be a lower limit due 
to the possible break-out of the SNR along the southern edge 
away from the molecular cloud to the north and east (Rho & 
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Fig. 2. Upper Left: NANTEN '-CO(J=1-0) image of thie W 28 region (linear scale in K km s"') for Vlsr=0 to 10 km s"' with VHE y-ray 
significance contours overlaid (green) — levels 4,5, 6cr as in Fig. 1. The radio boundary of W 28, The 68% and 95% location contours of 
GRO J1801-2320 and the location of the HIT region W 28A2 (white stars) are indicated. Upper Right: NANTEN ''CO(J=1-0) image for 
Vlsr=10 to 20 km s"' (linear scale and same maxima as for upper left panel). Bottom panels: Distribution of CO emission over the Galactic 
longitude and Vlsr plane integrated over Galactic latitude b ranges -0. 125° to -0.5° (left) and -0. 125° to -0.7° (right). The latter range is used to 
show the effect of extending the latitude range to encompass component A of HESS J1800-240. The bold circle indicates the suggested systemic 
velocity (7 km s"') of W 28 from the HI studies of Velazquez et al. (2002). 



Borkowski 2002). The HII regions, W 28A2 and G6. 1-0.6 are 
prominent in the 8.28 |J.m image (Fig. 3 right panel) from the 
Midcourse Space Experiment (MSX), showing that a high con- 
centration of heated dust still surrounds these very young stellar 
objects. 



5. Discussion 

Our discovery of VHE y-ray emission associated with dense 
(« > 10^ cm"^) molecular clouds in the W 28 field adds to the list 
of such associations after the detection of diffuse y-ray emission 
from the Galactic Ridge (Aharonian et al. 2006c), the associa- 
tion of HESS J1834-087 with the old-age SNR W 41 (Lemiere 
et al. 2005, Albert et al. 2006) and VHE emission discovered 
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Fig. 3. Left: VLA 90cm radio image from Brogan et al. (2006) in Jy beam"' (rebinned by a factor 1.2 compared to the original). The VHE 
significance contours (green) from Fig. 1 are overlaid along with the HII regions (white stars) and the additional SNRs and SNR candidates 
(with yellow circles indicating their location and approximate dimensions) discussed in the text. Right: ROSAT PSPC image — 0.5 to 2.4 keV 
(smoothed counts per bin from Rho & Borkowski et al. 2002). Overlaid are contours (cyan — 10 linear levels up to SxlO""* W m"- sr"') from the 
MSX 8.28 [im image. Other contours and objects are as for the left panel. The X-ray Ear representing a peak at the northeastern edge is indicated. 



from IC 443 (Albert et al. 2007). The VHE/molecular cloud as- 
sociation could indicate a hadronic origin for the parent multi- 
TeV particles where the y-ray emission (multi-GeV to TeV en- 
ergies) arises from the decay of neutral pions resulting from the 
interaction of accelerated protons (and higher Z nuclei) with am- 
bient matter of density n. In this case the y-ray flux would scale 
with cloud mass or density, and the total energy in accelerated 
particles or CRs penetrating the cloud(s). We note that a perfect 
correlation between the VHE and molecular cloud morphologies 
is not expected due to complex time and energy-dependent prop- 
agation of CR to and within the cloud (see Gabici et al. 2006 for 
a discussion). Projection effects are also likely to be important 
for the examples discussed here since the VHE emission could 
have contributions from clouds at different velocities, not nec- 
essarily physically connected to one another. For example the 
relationship between HESS J1801-233 and the W 28/molecu- 
lar cloud interaction is not entirely clear due to the overlapping 
molecular cloud components at Vlsr >10 km s"'. 

One should also consider accelerated electrons as the source 
of y-ray emission, via inverse-Compton (IC) scattering of am- 
bient soft photon fields and/or non-thermal Bremsstrahlung 
from the interaction of electrons with dense ambient matter 
Maximum electron energies however may be considerably lower 
(factor ~10 or more than that of protons) due to synchrotron 
cooling in magnetic fields and low shock speeds, in the ab- 
sence of strong electron replenishment. An assessment of the 
role of accelerated electrons requires consideration of the non- 
thermal radio and X-ray emission (where a convincing mea- 
surement of the latter is so far lacking), and also magnetic 
fields in this region. Such observations will also provide con- 
straints on synchrotron emission expected from secondary elec- 
trons resulting from primary hadron interactions with ambient 
matter (as discussed above). Relatively high magnetic fields 
B ~ 100(n/10'*cm"^)"^ [iG are inferred in dense molecular 
clouds (Crutcher et al. 1999). In addition, higher values are in- 



Table 2. Details for the molecular clouds towards the VHE sources in 
the W 28 field, assuming a distance d. 



VHE Source 




Vlsr 
(kms-') 


d 

(kpc) 






^^CR 


HESS J1801- 


-233 


0-25 


2.0 


0.5 


1.4 


13 


HESS J1801- 


-233 


0-12 


2.0 


0.2 


2.3 


32 


HESS J1801- 


-233 


13-25 


4.0 


1.1 


0.6 


23 


HESS JISOO- 


-240 


0-20 


2.0 


1.0 


1.0 


18 


HESS JISOO- 


-240 A 


12-20 


4.0 


1.0 


0.7 


28 


HESS J1800- 


-240B 


0-12 


2.0 


0.4 


2.3 


18 


HESS J1800- 


-240B 


12-20 


4.0 


1.5 


1.2 


19 



t Cloud mass xlO'' Mq 
t Cloud density xlO^ cm"^ 

§ Cosmic-ray density enhancement, ^:cr above the local value required to 
produce the E > I TeV VHE y-ray emission (using Eq.lO of Aharonian (1991)). 



dicated from Zeeman splitting measurements in the compact ar- 
eas (arcsecond scale) surrounding the 1720 MHz OH masers of 
the northeastern interaction region (Hoffman et al. 2005), co- 
inciding with HESS J1801-233. To the north of W 28, an- 
other potential source of particle acceleration is PSR J1801-23, 
where the VHE emission may arise in an asymmetric pulsar- 
wind-nebula (PWN) scenario (a primarily leptonic scenario), 
similar to HESS J1825-137 (Aharonian et al. 2006d). However 
with a spin-down power of E ~ 6.2 x 10^"^ erg s"' at distance 
d > 9.4 kpc, this pulsar appears unlikely to power any of the 
y-ray sources observed in the region. A PWN scenario would 
therefore require a so far undetected energetic pulsar 

In the case of a hadronic origin and following Eq.lO of 
Aharonian (1991), we can estimate the CR density enhancement 
factor kcR in units of the local CR density required to explain 
the VHE emission, given an estimate for the cloud masses and 
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assumptions on distance. Converting the VHE energy spectra in 
Table 1 to an integral value for E > I TeV, assuming distances 
of 2 and 4 kpc for the various cloud components, and that all 
the VHE emission in each source is associated with the cloud 
component under consideration, we arrive at values for A:cr in 
the range 13 to 32 (Table 2). 

Overall, these levels of CR enhancement factor would be ex- 
pected in the neighbourhood of CR accelerators such as SNRs. If 
the clouds were all at ~2 kpc, an obvious candidate for such par- 
ticle acceleration is the SNR W 28, the most prominent SNR in 
the region. Despite its old age, multi-TeV particle acceleration 
may still occur in W 28 (Yamazaki et al. 2006), with protons 
reaching energies of several lO's of TeV depending on various 
SNR shock parameters such as speed, size and ambient mat- 
ter density. In addition, CRs produced at earlier epochs have 
likely escaped and diffused throughout the region, a situation 
discussed at length in Aharonian & Atoyan (1996). Aharonian & 
Atoyan show for slow diffusion (diffusion coefficient at 10 GeV 
Dio ~10^^ cm^ s ' as might be expected in dense environments) 
CR enhancement factors in the required range could be found 
in the vicinity (within 30 pc - note that if at 2 kpc distance, 
HESS J1800-240 would lie ~10 pc from the southern circular 
boundary of W 28) of a canonical SNR as an impulsive accelera- 
tor up to ~ 10^ yr after the SN explosion (see their Fig. 1). In this 
sense, W 28 as a source of CRs in the region could be plausible 
scenario. 

The W 28 field however is a rich star formation region, and 
several additional/alternative sources of CR acceleration may 
be active. The SNR G6. 67-0.42 is positioned directly to the 
southeast of HESS J1801-233 (Fig. 3 left panel) while the SNR 
G7.06-0.12 is situated ~ 0.25° north of HESS J1801-233 and 
on the west side of the HII region M 20. M 20 itself may also 
be an energy source for the molecular clouds in this region. The 
SNR candidate G5.7 1-0.08 (Brogan et al. 2006) may also be 
responsible in some way for HESS J1800-24C given the good 
positional overlap between the two. These radio SNR/SNR can- 
didates are without a distance estimate making it unclear as 
to how they relate to the molecular clouds in the region. The 
morphology of HESS J 1800-240 displays several peaks, per- 
haps resulting from changes in cloud density and/or the pres- 
ence of additional particle accelerators and local conditions. For 
HESS J1800-24B, a potential energy source is the unusual ultra- 
compact HII region W 28A2 (G5. 89-0. 39), representing a mas- 
sive star in a very young phase of evolution. W 28A2 exhibits 
very energetic bipolar molecular outflows (Harvey & Forveille 
1988, Acord et al. 1997, Sollins et al. 2004) which may arise 
from the accretion of matter by the progenitor star The outflow 
ages are estimated at between ~ lO-' to 10** yr Recent observa- 
tions (Klaassen et al. 2006) suggest both outflows extend over a 
combined distance of ~ 2' (or ~ 1 .2 pc at c/ = 2 kpc), with total 
kinetic energy of 3.5 x 10"^^ erg. Surrounding the outflows is a 
very dense (> 10"* cm"-') molecular envelope of diameter 0.5' to 
1'. Despite the lack of any model to explain multi-TeV particle 
acceleration in such HII regions, its kinetic energy budget and 
its spatial overlap with a VHE source makes W 28A2 a tempting 
candidate for such acceleration. Already, there are two examples 
of VHE emission possibly related to the environments of hot, 
young stars — TeV J2032+4130 (Aharonian et al. 2005b) and 
HESS J1023-575 (Aharonian et al. 2007c). In this context, the 
HII regions G6. 1-0.6 and 6.225-0.569 may also play a simi- 
lar role in HESS J1800-24A. Among the prominent open clus- 
ters in the area, NGC 6523 and NGC 6530 ~ 0.5° southeast of 
HESS J1800-240, and NGC 6514 associated with M 20 ~ 0.7° 
north of HESS J 1801 -233 may also provide energy for CR pro- 
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Fig. 4. Energy fluxes of HESS J1801-233 and HESS J1800-240 (for 
regions defined in Tab. 1) compared to the £ > 100 MeV counterpart 
GRO J1801-2320. The power law fits and data points (summarised in 
Tab. 1) are also indicated: HESS J 1801-233 (solid blue line and points); 
HESS J1800-240 (open red points and solid line); GRO J1801-232 
(solid black points and grey Icr confidence band). 



duction. Finally, if the VHE emission is associated with truly 
distant cloud components approaching the Scutum-Crux arm at 
~4 kpc, undetected background particle accelerators would then 
play a role. 

Fig. 4 also compares the EGRET and VHE spectra. Given 
the degree-scale EGRET PSF GRO J 180 1-2320 remains unre- 
solved at scales of the VHE sources. Although the peak of the 
EGRET emission coincides with HESS J1801-233, we there- 
fore cannot rule out unresolved MeV/GeV components from 
HESS J 1800-240. Observations with GLAST wiU be required 
to determine the MeV/GeV components of the VHE sources. 



6. Conclusions 

In conclusion, our observations with the H.E.S.S. y-ray tele- 
scopes have revealed VHE y-ray sources in the field of 
W 28 which positionally coincide well with molecular clouds. 
HESS J1801-233 is seen toward the northeast boundary of 
W 28, while HESS J1800-240 situated just beyond the south- 
ern boundary of W 28 comprises three components. Our studies 
with NANTEN '^CO(J=1-0) data show molecular clouds span- 
ning a broad range in local standard of rest velocity Vlsr =5 to 
~20 km s"', encompassing the distance estimates for W 28 and 
various star formation sites in the region. If connected, and at a 
distance ~2 kpc, the clouds may be part of a larger parent cloud 
possibly disrupted by W 28 and/or additional objects related to 
the active star formation in the region. Cloud components up to 
~4 kpc distance (Vlsr >10 km s"') however, remain a possibil- 
ity. 

The VHE/molecular cloud association could indicate a 
hadronic origin for the VHE sources in the W 28 field. Under as- 
sumptions of connected cloud components at a common distance 
of 2 kpc, or, alternatively, separate cloud components at 2 and 
4 kpc, a hadronic origin for the VHE emission implies cosmic- 
ray densities ~10 to ~30 times the local value. W 28 could 
provide such densities in the case of slow diffusion. Additional 
and/or alternative particle accelerators such as HII regions rep- 
resenting very young stars, other SNRs/SNR candidates and/or 
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several open clusters in the region may also be contributors. 
Alternatively, if cloud components at Vlsr >10 km s"^ are at 
distances d ~ 4 kpc, as-yet undetected particle accelerators in 
the Scutum-Crux arm may be responsible. Detailed modeling 
(beyond the scope of this paper), and further multiwavelength 
observations of this region are highly recommended to assess 
further the relationship between the molecular gas and potential 
particle accelerators in this complex region, as well as the nature 
of the acclerated particles. In particular, further sub-mm obser- 
vations (eg. at high CO transitions) will provide more accurate 
cloud mass estimates, and allow to search for disrupted/shocked 
gas towards the southern VHE sources. Such studies will be 
valuable in determining whether or not W 28 and other ener- 
getic sources have disrupted molecular material at line velocities 
>10kms"^ 
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